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We discuss the results of a study of restricted solid-on-solid model hamiltonians for fee (110) 
surfaces. These models are simple modifications of the exactly solvable BCSOS model, and are 
Xf^ , able to describe a (2 x 1) missing-row reconstructed surface as well as an unreconstructed surface. 

■ They are studied in two different ways. The first is by mapping the problem onto a quantum 

0^ ' spin-1/2 one- dimensional hamiltonian of the Heisenberg type, with competing S^S^ couplings. The 

second is by standard two-dimensional Monte Carlo simulations. We find phase diagrams with the 
^ ' following features, which we believe to be quite generic: (i) two flat, ordered phases (unreconstructed 

D ' and missing-row reconstructed); a rough, disordered phase; an intermediate disordered flat phase, 

, characterized by monoatomic steps, whose physics is shown to be akin to that of a dimer spin 

state, (ii) a transition line from the (2 x 1) reconstructed phase to the disordered flat phase showing 
exponents which appear to be close, within our numerical accuracy, to the 2D-Ising universality 
class, (iii) a critical (preroughening) line with variable exponents, separating the unreconstructed 
phase from the disordered fiat phase. Possible signatures and order parameters of the disordered 
^ , flat phase are investigated. 
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I. INTRODUCTION 



On ■ Surfaces of fee metals, in particular (110) faces, display a variety of phase transitions which have been the subject 
of considerable experimental and theoretical work. Experimentally, the (110) faces of some fee metals - such as Au 
^ , or Pt - reconstruct at low temperature into a (2 x 1) missing-row (MR) or related structures, whereas other metals 
- like Ag, Ni, Cu, Rh, and Pd - retain (at least when clean) their bulk-like periodicity. As temperature is raised, 
' , \ reconstructed surfaces tend to show two separate transitions: a critical deconstruction transition, and, at a higher 
temperature, a Kosterlitz-Thouless roughening transition. OB] On the other hand, unreconstructed surfaces have not 



Q ' been shown, so far, to reveal a similar two-transition scenario. Only a roughening transition is well documented in 
O ' this case. ||,|| 

^ Based on theoretical considerations and on simulation work, an interesting and nontrivial interplay has been an- 

• I— I . ticipated between in-plane disordering, related to deconstruction, and vertical disordering, related to roughening, ||] 
' and many other studies have been devoted to the problem. P-p^ The situation is, in principle, somewhat different 
^ ] for the two types of situations, i.e., the unreconstructed and the MR reconstructed cases. On an fee (110) surface 
- - ' one can identify two interpenetrating rectangular sublattices, with origin, say, at (the "white" sublattice) and A at 
(-\/2x + y + z)a/2 (the "black" sublattice) where a is the lattice parameter, x = (001), y = (HO), and z — (110). The 
unreconstructed (ideal) surface has therefore two T — ground states, differing for the sublattice which occupies the 
top layer (see Fig. den Nijs has argued that, in such a case, the phase diagram should be qualitatively the same 
as that of a simple cubic (100) surface. |l^ In particular, den Nijs, |ll] Jug and Tosatti, Kohanoff et al., Q and 



Mazzeo et al. 11 1 argued that (110) surfaces like those of Ag and Pd(which do not reconstruct) are good candidates 
for realizations of preroughening, a critical (non- universal) transition from a low temperature ordered phase to an 
intermediate disordered fiat phase, previously identified in the context of restricted solid-on-solid models for simple 
cubic (100) surfaces. [|l^ In terms of the two ground states of the unreconstructed surface, the preroughening transi- 
tion can be viewed as a disordering of the surface due to the proliferation of monoatomic steps (see Fig. ||) separating 
terraces with one type of ground state from others where the other ground state is locally present. These steps retain, 
however, an up-down long range order - stabilized by a combination of up-up (down-down) step repulsion and entropy 
- which prevents the surface from being rough. [l5|] 

On the (2 x 1) MR reconstructed surfaces, with half of the (110) close-packed rows missing, the periodicity in the 
(001) direction is doubled. The surface has therefore four degenerate T = ground states, which can be classified by 
a clock variable 6* = 0, 7r/2, tt, 37r/2, according to the "color" and the position of the MR in the doubled unit cell, i.e., 
determined by which of the four sublattices sits in the top layer (see Fig. |[). The elementary extended defects which 
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one can consider here were discussed by Vilfan and Villain ||] and den Nijs |^ (see Fig. g). These are a) steps, which 
simultaneously change the average height by Ah = ±1, and the reconstruction variable 9 by A0 = tt/2 [clockwise or 
(3 X 1) steps] or A9 = — 7r/2 [anticlockwise or (1 x 1) steps], b) Ising wall defects with Ah = and A9 = tt, which 
can be seen as a tightly bound state between two steps of opposite sign (up and down), but same A9. ^ 

In this framework, den Nijs introduced a phenomenological four state clock-step model to describe the interplay 
between reconstruction and roughening degrees of freedom. |^ The model is formulated on a length scale larger than 
microscopic, through the introduction of a coarse grained lattice of cells in which a integer variable hr, representing 
the average height in the cell, and a clock reconstruction variable 9^ are defined. A bond in the lattice can be either 
empty (no defect), or occupied by an up or down step of either kind, or doubly occupied by an up and down step of 
the same kind (equivalent to an Ising wall), den Nijs found that when (1 x 1) and (3 x 1) steps have the same energy 

- the so-called zero chirality limit - the model displays two possible scenarios: (i) If the energy of an Ising wall is 
less than roughly twice the energy of a step Eg , temperature drives the system from the ordered phase to a disordered 
flat phase through an Ising transition, and then to a rough phase through a Kosterlitz-Thouless (KT) transition; 
(ii) When steps are energetically favored, > 2Es, the system undergoes a single roughening-plus-deconstruction 
transition, which is Ising-like for the reconstruction degrees of freedom and KT-like for the height degrees of freedom. 
The disordered flat phase present for E^^ < 2Es is quite clearly characterized by the proliferation of Ising wall defects 
(their free energy per unit length goes to zero at the deconstruction) . Accordingly, the surface shows a prevalence of 
9 = and 9 = tt terraces, say, over 9 = Tr/2 and 9 ~ in/ 2 ones. Using the terminology introduced in Ref. such 
a phase could be called disordered even flat (DEF). It has an obvious non-zero order parameter which counts the 
difference in the abundance of = 0, tt terraces over that of 6* = 7r/2, 37r/2 ones, and vanishing only in the rough phase. 
By contrast, when single steps dominate - i.e., 2Es < Ey^ - there is apparently no mechanism, in this simple model, 
which may stabilize the up-down long range order for steps, typical of disordered flat phases. It has been argued 
that suitable interactions penalizing the crossing of two up-up or two down-down steps - not considered by den Nijs 

- could stabilize such a hypothetical step-dominated disordered flat phase. A disordered flat phase of this kind 

- termed DOF in Ref. [|l2[ - should be characterized by an equal abundance of all types of cells 9 = 0, 7r/2, tt, 37r/2, 
i.e., by a vanishing of the order parameter characterizing the DEF phase. [Il2| 

Interestingly, the situation does not change much in the so-called strong chirality limit considered by den Nijs, i.e., 
when anti-clockwise steps, say, are very costly and thus completely negligible. In such a case the problem may be 
mapped onto a one-dimensional fermionic model containing a Hubbard type on-site step-step interaction U such that 
the energy of a Ising wall configuration (doubly occupied site) is Eyj = 2Es + U . For [/ < 0, the result is the same 
as in the zero chirality limit case, i.e., a DEF phase is obtained. When steps dominate - i.e., for [/ > - two distinct 
rough phases appear, the deconstruction transition is no longer of the Ising type, but, again, no disordered flat phase 
exists. 

The variety of possibilities introduced by the phenomenological models is thus very rich. Do microscopic models 
display just the same, or any new features, one might ask? In the light of the previous discussion, the question 
naturally arises as to what kind of disordered flat phase (or phases) is realized in simple solid-on-solid (SOS) models. 
The question has been considered by Mazzeo et ai, who have introduced, and studied by Monte Carlo simulation, 
a restricted SOS model which is able to deal with both unreconstructed and reconstructed situations. ||ll|] Their 
model - which we will refer to as "i^a-model" - is a simple modification of the exactly solvable body-centered 
solid-on-solid model (BCSOS) obtained by adding a further neighbor interaction, which can stabilize the (2 x 1) 
MR reconstruction if required. For a reconstructed case, they find two-transitions: a 2D-Ising deconstruction to a 
disordered flat phase, and a KT roughening at a slightly higher temperature. The unreconstructed case studied also 
shows two transitions, with a non-universal critical transition to a disordered flat phase followed by a KT roughening. 
[ pl| The nature of the disordered flat phase present in the model was, however, not fully characterized. 

In the spirit of investigating simple but fully microscopic models, Santoro and Fabrizio have studied a slightly 
different modification of the BCSOS model, which will be referred to as "X4-model" . ||l^ They found that this model 
could be mapped onto a quantum spin-1/2 Heisenberg chain with further-neighbor interactions. (l9j The phase diagram 
they obtained has two low-temperature ordered phases ~ unreconstructed or (2 x 1) MR reconstructed, depending on 
the parameters of the interactions - a high-temperature rough phase, and an intermediate disordered flat phase. The 
physics of the disordered flat phase was argued, by analytical arguments, to be akin to that of the dimer quantum 
spin phase studied by Haldane, pO| i.e., a doubly degenerate state which breaks translational invariance and in which 
dimer-dimer correlation functions acquire long range order. 

In this paper we extend and apply the approach of Ref. ]l8[ ] in such a way as to provide a unified picture of the 
phase diagram of all the simple BCSOS-like microscopic models of fee (110) surfaces studied so far. First we show 
that a quantum spin-1/2 hamiltonian also underlies the more general restricted SOS model where both the couplings 
considered in Refs. and are included. The spin-1/2 model is, in all cases, a Heisenberg chain with 5f5| 
competing antiferromagnetic interactions ranging up to third neighbors. Secondly, we unambiguously show that the 
dimer-phase scenario is realized in the disordered flat phase, and discuss in detail the surface physics implications 
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of such a scenario. Long-range dimer order suggests (see section ^ a particular type of long-range order for the 
correlation function between (2 x l)-steps, and also between local surface maxima. In particular, one is led to study 
an order parameter Pbw - previously introduced by Bernasconi and Tosatti ||l2| - which measures the difference in the 
number of local surface maxima belonging to the white and the black sublattice of a fcc(llO). Due to "shadowing" , or 
to surface lattice contraction, this order parameter is related to antiphase scattering intensity of He-atoms of X-rays, 
respectively, and is thus a quantity of direct interest. (See section 0.) Pbw is studied by finite-size analysis of 
classical Monte Carlo data, and found to be non-zero in the disordered flat phase of both the K3- and Ki-model. It 
has a non-monotonic temperature behavior, vanishing only at preroughening and in the rough phase. 

The present paper is organized as follows. Section |l| introduces the BCSOS-type of models which we consider. In 
Section III we show in detail how these models may be mapped onto quantum spin-1/2 chain problems. In Sections 
IV and ^wc discuss the phase diagrams as well as the physics of a "dimer" disordered flat phase. In Section VI 
we present the results of our Monte Carlo simulations and discuss possible experimental signatures of a dimer phase. 
Section VII, finally, contains some concluding remarks as well as a discussion of open problems. 



II. RESTRICTED SOLID-ON-SOLID MODELS FOR A FCC (110) METAL SURFACE 



The (110) surface of a fee lattice is comprised of two interpenetrating rectangular sublattices of lattice constants 
Qx = V^Uy, which we will conventionally refer to as the white (W) and the black (B) sublattice. The surface lattice 
basis vectors are x = Oa-x and y = Uy-y, where x — (001), and y — (110). In the ideal unreconstructed (110) surface, 
one of the two sublattices lies above the other at a distance — ay/2. Within a solid-on-solid (SOS) framework, ]2l| ] 
one associates to each site r of the lattice a height variable which can take only integer values (take = 1). The 
models we are going to study have an additional restriction, in that the height difference between each site and its four 
nearest neighbors (belonging to the other sublattice) are forced to be Ah = ±1. A height difference of is therefore 
excluded, as well as larger values of Ah (they are energetically more costly). As a consequence, the values of hr are 
forced to have opposite parity on the two sublattices, say even on W and odd on B. This restriction is typical of the 
BCSOS model. [|l^ It is probably justified for a metal, where strong inward relaxation makes the first-second-layer 
bonds extra strong. On the contrary, it should not be expected to hold for, say, a rare gas solid (110) face, where 
these bonds are in fact weaker. 

Our hamiltonian is written as 



Ti. = "Hbcsos + AH , (1) 

where Hbcsos describes interactions between sublattice nearest neighbors, and ATI takes into account further- 
neighbor interactions which will favor or disfavor reconstructed phases. Specifically, Hbcsos is given by 

HbCSOS = K2y J2{hr+y - Kf + K2^ ^(/ir+x - Kf , (2) 
r r 

with different coupling strengths in the two directions to take into account the anisotropy of the surface. K2y will be 
always assumed to be positive and is generally the largest energy in the problem. The correspondent physics is that it 
is very costly to create a kink on a tightly packed (110) row. The absolute value of K2x, i-e., of the coupling between 
rows, is instead much smaller, since atoms in neighbouring rows are only second bulk neighbours. For K2x > 0, the 
(110) surface is stable in its (1 x 1) unreconstructed form. If AH = we recover the BCSOS model which is exactly 
solved through a mapping to the six vertex model, [ pT| and shows a single transition. This is of the Kosterlitz-Thouless 
(KT) type, between a low-temperature ordered (unreconstructed) flat phase and a high-temperature disordered rough 
phase. For K2x < the surface becomes unstable against (110) step formation. In this case the final stable state is 
determined by more distant interactions, contained in AH. As for AH, two possible simple choices have been made 
in the literature, corresponding to what we will refer to as the "i^s-model" and the "Ki-modeV . The i^a-model has 
been introduced by Mazzeo et at, pl and is defined by 

A'H(if3) ^ X! [(/ir+x+b - hr)^ + (/ir+x-y+b " Kf] , (3) 

r 

with K3 > 0, and b = (x -f y)/2 (see Fig. |^). The introduction of this further-neighbor interaction stabilizes the 
(2 X 1) MR reconstructed phase. ||ll|] In fact, it is very easy to check that K2x < 0, > stabilizes an ordered 
succession of up and down (110) steps, which is precisely the (2 x 1) MR state. An alternative way of stabilizing the 
same (2 x 1) MR state against (111) faceting is realized with the K^-model, whose AH reads 
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r 



(4) 



with Ki > 0. The fourth neighbor interaction in the x-direction has the effect, once again, of increa sing the energy of 
configurations with |/ir+2x — h^] =4. This model was originally proposed by Kohanoff and Tosatti, |2^] and has been 
recently studied in detail in Ref. More generally, we could include both types of couplings by taking 

An = AH(K3) + AH(K4) • (5) 

In subsequent calculations and simulations, the lattice will be taken to have N,. = Nj; x iV^ primitive cells, i.e., 
2Nx X Ny sites. Periodic boundary conditions are assumed in both directions. A schematic representation of the 
lattice and of the interactions considered is given in Fig. ||. 

The classical T = ground states for both models are easy to work out as a function of the dimensionless ratio 
IC = K2X/K3 or /C = K2x/Ka- For both models one finds that /C > corresponds to an unreconstructed ground state 
[or (1 X 1)], whereas for —4 < /C < the ground state is (2 x 1) MR reconstructed. For the sake of completeness, 
we mention that, for the i^a-model, the ground state degenerates into an infinite (111) large facet as soon as /C < —4. 
PI For the iir4-model, on the contrary, an infinite (111) facet sets is only for K. < —8, and there is a whole series of 
intermediate regions [(12 — 8n)/{n — 1) < K. < (20 — 8n)/{n — 2) with n > 3] where the ground state is {n x 1) MR 
reconstructed. In the following we will be interested exclusively in the region of parameter space where the interplay 
between unreconstructed and (2 x 1) MR reconstructed phases takes place, hence K, > —4. 

Neglecting adatoms, vacancies and (as necessary in SOS models) overhangs, the defects which should play a role 
in the disordering and roughening transitions are unbound steps and bound pairs of steps, i.e., Ising domain walls. 
Fig. ^ illustrates the most relevant defects for both a (1 x 1) and a (2 x 1) surface. The ground state energies of 
these defects are given, for both the K3 and the iir4-model, in the table. It is worth noting that the ifs-model has 
defects whose energy goes to zero as K2X 0. These are the (2 x l)-step and the (2 x l)*-wall in the (1 x 1) 
case, and the anticlockwise-step and the Ising wall in the (2 x 1) reconstructed case. Since this leads to zero-point 
entropy, the K2X — > region is therefore somewhat unphysical for the iCg-model, where one might expect disorder 
to occur at very low temperatures. [ p3| As a second point, we observe that in the if4-model a combination of two 
Ising walls wins against combinations involving clockwise-steps for /C — K2x/Kn > —1, while it always wins against 
anticlockwise-steps. In principle, therefore, a DEF (wall dominated) phase seems to be plausible for — 1 < /C < in 
the i(r4-model. Later on we will present results which show how ground state defect energy considerations can be 
somewhat misleading: the disordered phase obtained does not have the features of an ideal DEF. 



III. MAPPING INTO A QUANTUM SPIN-1/2 CHAIN 

An elegant and convenient way of studying the temperature phase diagram of the classical two-dimensional models 
introduced in the previous section consists in mapping them into one-dimensional quantum problems, [ p^ The general 
procedure is well known, [ p4[jl9|j25[ | but we review it here for our specific case, and for the reader's convenience. 

The method consists in viewing the y-direction of the lattice as the (imaginary) time direction of an appropriate 
ID quantum problem, different row configurations in the x-direction being viewed as subsequent time-slices for the 
quantum problem. The physical requirement which will turn out to be important is that the coupling K2y in the 
y-direction is strong, while the other couplings are much weaker (strong anisotropy limit) . The starting point for the 

mapping is a transfer matrix formulation of the classical partition function Z. The notation r = (i, j) (and hr = h!f^) 
for the lattice sites used in the present section is illustrated in Fig. |^: the value of j, the time-slice index, is the same 
along each dashed zig-zag line shown in Fig. l]; within each zig-zag line, the W and the B sublattices are characterized, 
respectively, by even and odd values of i. The classical partition function is given by 

where \h^^'^) — {h!f^ : i — 1, ■ ■ ■ ^ N^} is the j-th row configuration (a dashed zig-zag line in Fig. |^, containing Ng — 2Nx 
sites), and T is the classical transfer matrix. Periodic boundary conditions have been used in the y-direction, and 
are understood in the x-direction. It is also understood that configurations differing by a uniform shift of the heights 
should be included only once in the partition sum. For the models we are considering, the transfer matrix elements 
read: 

^ B„„ exp{-pK2yJ2{h^'''^-h^^r -l3K2xY.{h^2-l^^f}, (7) 
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where i?„„ is the Boltzmann weight contribution due to further neighbor interactions. For the ifs-model, _B„„ is 
given by: 



i even 

X exp{-^ f; - h!f^f + if^^ - /.pV]} , (8) 

i odd 

whereas for the i4r4-niodel the result is: 

= exp{-/3X4 V(/.!i^4 - f^^n ■ (9) 



For the model in which both couplings are present, one clearly has, 

i?„„ = Bi^^^i^i^^) . (10) 

Notice that in the partition function, Eq. ^, it is implicitly assumed that the configurations included have to fulfil the 
BCSOS constraint Ah = ±1 for nearest neighbors. As a consequence, within each row we must have h['2i — h\^'^ = ±1. 
Therefore, we can associate to any row configuration \h'^^^) a state \j) in the Hilbert space of a quantum spin-1/2 
chain (of length Ns) by the relationship 

S. = \{h^,~h^^)- (11) 

(In doing so we lose information on the absolute height of the surface, which is however irrelevant in a static context.) 
Fig. H illustrates the explicit mapping of the (1x1) and (2x1) ground states in terms of spin configurations. 

The idea is now to reproduce the Boltzmann factors appearing in the matrix elements of the classical transfer 
matrix {h'^^^^^\T\h'^^^) by a suitable quantum operator Tq in the spin Hilbert space, i.e., 

= (.? + 1|Tq|j) , (12) 

where \j) and |j + 1) are the quantum spin states corresponding to \h'^^^) and \h^^^^^), respectively. In certain cases, 
the exact expression for the quantum operator Tq can be worked out quite easily. Tq for the iir4-model has been 
derived in Ref. ||l8|| . The exact Tq is, however, of no practical use, being a product of non-commuting terms involving 
spin-1/2 operators. The crucial step which makes the whole mapping useful is the so-called time- continuum limit 
or strong anisotropy limit. Physically, one assumes that the "time" direction coincides with the "hard" direction of 
the classical problem, i.e., that the coupling in the y direction is much stronger than the other couplings. This is 
plausibly the case for the fcc(llO) surface, where the y = (110) direction is hard, and the x = (001) direction is 
soft, as discussed above. Anisotropy, moreover, is not expected to play any role in the qualitative shape of the phase 
diagram. |l^ 

In the strong anisotropy limit the quantum operator Tq will reduce to the imaginary-time evolution operator e~'^^^ 
for a suitable hamiltonian Hs, with r ^ 0. To find Hs, assume /9i^2y to be large, so that e~'^^^^^ ~ (J/2)t, t being 
a small quantity (and J of order one). Assume also all the other couplings to be small and proportional to r, i.e., 
(3K2x oc r, (3K^ oc r, and (3K4^ cx r. We need a spin hamiltonian Hs such that Eq. ^ is verified with 

TQ^e--"' ^1~tHs+0{t^) , (13) 

up to first order in the small quantity r. A diagonal matrix elements of T reads, using Eqs. to 0, 

(/.(^■)|f|/.(^)> = l-pK,xY.^2-h^^f-PK^(h^,-h^^f 

i i 

-PK,J2{h!^,^h^^r +0{r') (14) 
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where we have expanded all the exponentials up to first order in small quantities (oc r). The first requirement for Hs 
is that its diagonal terms must give the same result, i.e., using Eqs. |lj and ^ 

{j\Tq\j) « 1^t{j\Hs\j) + 0{t') 

i i 

-A(3KiJ2^S, + + S,+2 + ^.+3)' + O(t') . (15) 

i 

An off-diagonal matrix element of T must contain a Boltzniann factor e~'^^^^y = (J/2)t, for each site i such that 
/ip^^-* = /ip'' ± 2. Therefore, up to first order in r we need to consider only row configurations which differ 

from h^^') only at a single site i. Let {• • • , hi-i, hi, /i^+i, • • •} be the local configuration of row j around such a site i. 

It is easy to realize that, in order to satisfy the BCSOS constraint with h^^^^ = h'f^ ± 2, and /i^"'^^'' = h^j^^ for k i, 

(?) (1) (1) 
the only possibility is to have hY_i ~ hi^i = ± 1, i.e., 

\h^^^) ^{■■■,h,±l,h,,h,±l,---} ^ \h^^+^^) ={■■■, h,±l,h,±2,h,±l,-- ■} . 
In terms of the corresponding spin configurations, this simply leads to a spin flip of the spins at sites i — 1 and i, 

|j) - I • • • , = T^, 5. - • • •) ^ b + 1) = I • • • , = ±i, S., = T^, • • •) • 

The corresponding T-matrix element reads, up to first order in r, 

«-rO- + l|i/5|i) . (16) 

It is easy to verify, in conclusion, that the correct form of Hs verifying Eqs. |l^ and |l6| is given, neglecting an overall 
constant, by 

2 = 1 1 = 1 

where the spin couplings are related to the original couplings as follows: 

Tj = 2exp (-4/3X2^) 
r J, = mK2^ + 2K3 + 3Ki) 
TJ2 = 8p{K3 + 2Ki) 

TJ3 = 8PK4, . (18) 

Indeed, the spin-fiip part of Hs reproduces the off-diagonal matrix element in Eq. whereas the terms give 

rise to the correct diagonal matrix element in Eq. 

It is well known that this kind of mapping is such that the free energy per site of the classical problem - given by 
the maximum eigenvalue of the transfer matrix - is related to the ground state energy per site of the one-dimensional 
quantum problem, (3f — recs- The temperature clearly enters through the spin couplings, see Eq. |l^, so that any 
genuine singularity of the classical free energy versus temperature can be seen as a ground state energy singularity 
for the quantum problem as a function of the couplings Jz/J, J2/J and Js/J- Moreover, temperature averages for 
correlation functions of the classical problem can be likewise rewritten in the form of ground state averages for the 
corresponding quantum correlation function. p4| ] In summary, to obtain information about the temperature phase 
diagram of the classical model one studies the ground state phase diagram of the spin chain model, 

Before entering into the discussion of the phase diagram, let us clarify that the quantum mapping not only gives 
the correct critical behavior of the transitions (if anisotropy is not "relevant" in the rcnormalization group sense), 
but provides also quantitative results on the transition temperatures that are expected to be quite good even if the 
anisotropy is in reality only weak. As a simple check of this point, consider the exactly solvable BCSOS case, whose 
transition temperature is given by |]l7| 
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The BCSOS model maps - see Eq. 18 



KT transition at the isotropic point Jz 
pi'^'' satisfying 



into the nearest neighbor XXZ Heisenberg chain, which is known to have a 
= J. po|] Using Eq. |l^, this readily implies a predicted transition temperature 



(20) 



Fig. ^ shows both the exact (solid line) and the quantum mapping transition temperature (dashed line) for the 
BCSOS as a function of the anisotropy ratio K2x/K2y The results agree within a few percents even for rather weak 
anisotropics, such as K2x/K2y = 0.2, and remain reasonable all the way to the full isotropic case, K2x/K2y = 1. 



IV. PHASE DIAGRAMS 



A. Phase diagram of the spin-1/2 chain 

The spin chain hamiltonian corresponding to our modified BCSOS model is a Heisenberg XXZ model with a second 
neighbor and a (less important) third neighbor S^S^ coupling 

J N N 

i=l 1=1 

Qualitatively - and for not too large values of the couplings Jz, J2, J3 ^ the physics of such a model is closely related 
to that of the Heisenberg chain with spin-isotropic second neighbor interactions 

-ff = - + Sr-S.^i] + "^yzS^'S-^i + J2Si ■ Si+2\ , (22) 

i i 

which we will refer to, in the following, as the J — J2-model. ||2^ Haldane has discussed the qualitative phase diagram 
of the J — J2-model in the context of the Luttinger liquid framework, with special emphasis on the role played by 
umklapp processes in the underlying spinless fermion problem. pO| (For a detailed quantitative analysis see Ref. [^) 
For the spin model in Eq. |2^, the discussion goes along similar lines. For completeness we will give in subsection pV C| 
some details of this analysis based on standard techniques of one dimensional systems. 29 



Even a simple mean field treatment, however, is quite instructive about the nature of the ordered phases which 
are to be expected. The starting point is to perform a Wigner- Jordan transformation from spin variables to spinless 
fermion operators i.e., = c\ci — 1/2, 5*^^ = cj exp (*7r^^^^ rij). Neglecting constants and terms proportional to 
the total number of fermions, the spin model is then rewritten as the following spinless fermion model 

N N 

Hp = ^[4c.i+i + c|_|_iCj] + ^[J^riiTii+i + J2nini+2 + J3"-i"-i+3] 

t=l i=l 
BZ ^ BZ 

= ^<^kclck + — ^v{q)piq)p{-q) , (23) 

k q 

where Ck — N^^^^ e^^^°'^Cj, with k belonging to the first Brillouin zone [— 7r/a,7r/a] {BZ), and p{q) is the Fourier 

transform of the fermion density operator, p{q) = J2k'^\'^k+q- Here ek ~ —2tcos{ka) with t — J/2, and v{q) is 
the Fourier transform of the interaction potential, v{q) = Jz cos (qa) + J2Cos{2qa) + J3Cos{3qa). Since J2i^i — 
rii — N/2, zero total magnetization for the spin system implies half-filling for the fermions, i.e., a Fermi surface 
consisting of two Fermi points at ikp, with kpa = tt/2. In absence of interaction (Jz = J2 = J3 = 0, i.e., the XY spin 
chain) we have a simple free-fermion problem. The two Fermi points induce nesting with a wave- vector 2kp ~ ir/a, 
a hint that the system would tend to open-up a gap at the Fermi surface by developing long-range order (LRO) with 
wave-vector ir/a and making the average Ak = different from zero. A standard mean-field factorization of 

the quartic term in Eq. assuming (c^_|_^y^Cfe) 7^ 0, leads to the following mean-field hamiltonian (we take a — 1 
from now on) 

RBZ 

Hf^ = h[cick - cl+,cfc+^] -I- [Afe4cfe+, -I- H.c] , (24) 

k 
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where h = Ck - Efc'esz "(^ " k'){clck'), and Afe = {2/N)J2k'eBzH^) " " k')]{cl+,Ck'). Here RBZ 

stands for the reduced magnetic Brillouin zone (— 7r/2, 7r/2). Diagonahzing the simple 2x2 problem appearing in 

, one readily finds that the "anomalous" average A^. = (c^^^cj.) is simply related to A^, — —Ak/{2Ek) with 

Ek = \/(\^ l^fcPi ^^^d that efc and A^ have to obey the self-consistency conditions 

^ RBZ 

e, = ek + -J2 [^(^ - - - ^' + '^)]f7 
fc' 

^ RBZ 

= "TV E ;^{[^(^) - - ^')]^'^' + - - ^' + • (25) 

Let us consider, for definiteness, the case J3 — 0. Solving the self-consistency equations one finds that if v{t:) — 
J2 — Jz is sufficiently negative (J2 less than « 0.4J^), A^ is real and has the form A^ = (5o + (52COs(2fc). Upon 
increasing J2, a transition occurs to a phase in which A^ is purely imaginary, A^ = i5i sin (fc). The transition appears 
to be first order in mean-field. To understand the meaning of the two phases, consider the average values of the 
fermion density (ri^), and of the bond kinetic energy (c|ci+i -I- c|^-^Ci). A simple calculation shows that 

RBZ „ . 
RcAfc 

N ^ Ek 

2 ImA 

(clc,+i + 4^,c,} = const + (-1)' - ^ sin (fc)-^ . (26) 

The phase with A^ real (small J2) is therefore a 2kp site-centered charge density wave (CDW) (i.e., a Neel phase with 
titi LRO in the spin language, or an ordered fcc(llO) face in the surface language), whereas the phase at larger J2, 
with Afc purely imaginary, is a 2fci;' bond-centered charge density wave. We can picture the latter by assuming that 
on every other bond the fermions are in a state which maximizes the kinetic energy of the bond, 

*Bond-CDW = n , (27) 

i even 

as opposed to the ideal site-CDW state (the Neel state titi) 

*Site-CDW= n 4|0) • (28) 
i even 

The spin interpretation of the bond-CDW state is quite obviously a dimerized spin state with every other bond engaged 
in a singlet, (ti — it)- ll3 Unhke the site-CDW, where every second neighbor is occupied and pays an energy J2, a 
bond-CDW reduces the second-neighbor average occupancy to about 1/2, and is thus favoured upon increasing J2. 
As will be discussed in detail, this spin dimcr phase corresponds to a disordered flat phase in the surface language. 
Clearly, for very large J2, the system will eventually prefer to minimize second-neighbor occupancies by forming a 
site CDW of double periodicity (i.e., a fc_F-CDW), which we can picture as 

*kp-CDW= n 44i|0) ■ (29) 

i—4n 

(Such a state corresponds to TTii LRO in spin language, or a (2 x 1) MR reconstructed face in the surface language.) 
This phase can be included in a mean-field treatment by allowing, in the factorization of the quartic term, anomalous 
averages of the type (c^^^yjCfe); as well as the previous one (c[._(^^Cfe). 

A finite-size scaling study of the spin model readily confirms most of the qualitative features of the mean-field 
phase-diagram. A quantitative phase diagram for the spin model corresponding to the iCa-model, i.e., Eq. ^ with 
J3 = 0, is presented in Fig. |^. The procedure to obtain such a phase diagram from a finite-size scaling study of chains 
up to iV = 28 sites, pQ] was described in detail in Ref. flS]. (See also Ref. PSI.) A similar phase diagram for the spin 



chain corresponding to the if4-model, i.e., Eq. ^with J3 = J2/2, was presented in Fig. 1 of Ref. |18|. For the purpose 
of a general discussion, we reproduce in Fig. ^ the essential qualitative features of the spin chain phase diagram for a 
generic J3 — aJ2 with < a < 1/2. The model has a spin liquid XY-like phase at small J2 and J^, which corresponds. 
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in the fermion language, to a spinless Luttinger liquid characterized by a certain Luttinger exponent K. (See IV C for 
more details on this discussion.) At a given universal value of the Luttinger exponent {K = 1/2), the Luttinger liquid 
phase becomes unstable - because of umklapp processes - against two different (gapped) ordered phases, depending 
on the sign of the effective coupling of the umklapp term: a Neel phase with titi LRO, for large and small J2, and 
a dimer phase, for larger J2. Both phases have a gap in the excitation spectrum, and a doubly degenerate ground state 
which breaks translational symmetry. |20 2^ These two ordered phases are separated by a critical line of continuously 



varying exponent, labelled 'PM' in Fig. 3[along which the effective coupling of the umklapp term vanishes and the 
system behaves as a Luttinger model with a Luttinger exponent 1/8 < if < 1/2. Beyond the point 'M' in Fig. the 
nature of the line changes from non-universal to first order. For even larger values of J2 the other ordered phase, with 
spins acquiring ttii LRO and a fourfold degenerate ground state, wins over the dimer phase. This is the only feature 
of Figs. ^ and ^ which is qualitatively new with respect to the phase diagram of the J — J2-model (Eq. ^2|). pl| ] 

The nature of the line separating the ttii phase from the dimer phase is an open issue. Previous studies of the K^- 
model pl| and of the spin chain found exponents which appeared to be compatible with the 2D-Ising universality 
class. Recently, a transfer matrix study of a 2D model closely related to the isr4-model has found exponents which 
are incompatible with Ising. |^2| A definite answer to the nature of this transition, possibly connected to the presence 
or absence of the multicritical point M in the phase diagram, [ ^2| needs further study. In spite of this uncertainty, we 
will continue to refer to this line, for convenience, as "Ising" . 

A second open issue concerns the region of the phase diagram where the KT line and the "Ising" line seem to 
approach each other. A relevant question, which we have not been able to answer, is whether the KT and the "Ising" 
lines actually merge, and, if so, what is the nature of the resulting line. 



B. Phase diagram of the modified BCSOS models 

The translation of Fig. || into a temperature phase diagram for ii'a-model, using Eqs. |lj, is shown in Fig. |^. p^ ] 
The generic phase diagram of our modified BCSOS model, in the {T,K2x) plane for given values of K3 and K4, is 
qualitatively sketched in Fig. |8[ ||2^ Four phases are found in a region of parameters relevant to the unreconstructed 
and (2 x 1) MR reconstructed case. At very high temperatures, there is a rough phase. It corresponds, in the 
spin problem, to the region close to XY-modcl point (J^ = J2 — J3 = 0) in which spin-sp in co rrelation functions 
decay as power laws at large distances (the Luttinger liquid or Gaussian model, see subsecti on |IV C ) . A large-distance 
uniform term of the type —K/{2'K'^n'^) in the spin-spin correlation function {S§S^) - see Eq. [11|- implies a logarithmic 

divergent height-height correlation function G{n) = {[hl^^^ — /lo^"*]^), signalling a rough phase. Indeed, using Eq. |ll| 
and translational invariance of the spin-spin correlation function, one verifies that 

G{n) = ([/.(°) - 4")]^) 

n— 1 n 

= 4 {StS^) = n + 8 ^(n - i){S^,S^;,) 

i.j—0 i—1 

4K 

= — — In (n) + • • • , n — > 00 . (30) 

At low temperatures, corresponding to large values of Jz/J and/or J2/J in the spin-chain problem, a (1 x 1) and a 
(2 X 1) ordered phase are present for K2X > and < 0, respectively. The (1 x 1) and (2 x 1) ordered phases 
correspond, respectively, to titi and ttii LRO for the spins (see Fig. The other phase appearing in Figs. |^ 
and ^ is a disordered flat (DF) phase. It corresponds, in the spin language, to the dimer phase (see Section ^ for 
a more extensive discussion). The transition line from the (2x1) reconstructed phase to the DF phase is labelled 
as "Ising", in spite of the fact that its nature is not completely assessed (see previous section). The critical line 
separating the unreconstructed phase from the disordered flat phase has variable exponents: it is the preroughening 
line. The parameter K appearing in Eq. ^ is the Luttinger exponent. In the rough phase K > 1/2. Along the 
preroughening line correlation functions still behave as power laws with exponents related to K; Eq. BO is still valid, 
with 1/8 < A' < 1/2. 



C. Spinless Luttinger liquid and the variable exponent line 

We now discuss in more detail how to extract, using standard techniques of one-dimensional systems, a qualitative 
phase diagram for our spin chain model and some useful information about the variable exponent line. The reader not 
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interested in technical details might jump directly to the next section, where the surface interpretation of the dimer 
spin phase is discussed. 

The starting point is the spinless fermion model in Eq. |2^. The low-energy physics of such a model, as long as the 
interactions are not too strong, can be conveniently analyzed by going to the continuum limit, a — > with Na = L 
fixed. One linearizes the fermionic band around the two Fermi points at ±kp, and introduces a right (p — +) and 
left {p = — ) branch of fcrmions, with fields "tppix). All the interactions processes in which particles are scattered 
in the vicinity of the Fermi points are then classified in the so-called "g-ology" scheme. The resulting continuum 
fermionic model TLp turns out to be a sum of two terms 

'Hp — ^Luttingcr + Wumklapp j (31) 

where T^Luttingor is a spinless Luttinger model, |Q 

"^Luttinger ^'"F^^ l dx : i/'p (a:) [-ipV - kF]'4'pix) ■■ 

+ [54V,p +52V -p] / dx : pp{x)pp'{x) : , (32) 

"T^j. Jo 

P,p'=± 

and Humkiapp represents the crucial umklapp processes, i.e., processes where two fermions are scattered from the 
vicinity of one Fermi point to the opposite one, 

^umklapp =53 / dx[: ^pl_{x)tp_{x)tpl_{x)ip_(x) : +H.C.] . (33) 



Jo 

(Umklapp processes would not conserve the momentum for a general filling: at half-filling, however, momentum 
conservation is fulfilled modulo a reciprocal lattice vector, G = Akp — 2tt.) Here vp = 2t = J is the Fermi 
velocity, and Pp{x) —: ^pj,{x)ipp{x) : is the density operator for the p-branch of fermions. (The : • • • : stands for a 
normal ordering procedure, as explained in [^ .) Neglecting lattice renormalization effects we have, for the Luttinger 
couplings (74 = f (0) = (J^ + J2 + J3) and 172 — ^{0) — vin) — 2(Jz -I- J3), whereas the umklapp coupling reads 
33 = vin) = {—Jz + J2 — Js)- We stress the important point that the sign of umklapp coupling g^ results from a 
competition of Jz and J2. We will see that this fact is crucial to the existence of a line with variable exponents. 

The final step is to hosonize the hamiltonian in Eq. This is achieved by introducing a bosonic representation of 
the fermionic fields p3| , p9[ 

iPp{x) = _^rype*f'=^^e*'"^''(^) , (34) 
v27ra 

where a is a short-distance cut-off, and rjp = rj^ are Majorana fermions (77^ = 1) ensuring correct anticommutation 
properties among right and left-moving fermions. The field 4>p{x) is related to the fermion density as follows 

Pp{x) = : ^ljl{x)il;p{x) : = -^Vcjjpix) , (35) 

and is expressed in terms of standard boson operators bp{q) as: 

1/2 



q>0 ^ 



(Here q = (27r/L)n, with n integer.) The continuum model in Eq. pT| translates, in bosonic variables, into a quantum 
sine- Gordon problem J53[|2^] 

(36) 

where we have introduced the canonical field $(a;) — {(f>+ + (f>-)/^/4:n and its conjugate momentum H(.t) = — V(^!)+ — 
(/)_)/'\/47r. The Luttinger model, Hmttingor, is equivalent to free bosons i.e., to a massless Klein-Gordon (or Gaussian) 
problem {V = 0). Q K is the crucial parameter governing the low-energy physics of the problem, vs being the 



HsG^Y J '^'^ + ]^(V*)'] + -^^^ j dx cos (Vl6^$) , 
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velocity of the sound-like gapless excitations. In weak-coupling we have vs/K = J[\ + (3 -I- J2 -I- 3J3)/ttJ -I- • • •] , 
and vsK = J[l — {Jz — J2 + J3)/ nJ+ ■ ■ ■]. The umklapp term in Eq. |3^, rewritten using Eq. ^ gives rise to the cosine 
term, with V — —2gj, = 2{Jz — J2 + J3) + ■ ■ ■■ The renormalization group equations for the sine-Gordon problem are 
well known, and have the Kosterhtz-Thouless form 

^^2[1^2K{l)]Vil), (37) 

with V simply proportional to 1^. V ^ with K > 1/2 is a line of stable fixed points which represent the Luttinger 
(or Gaussian) model. V = with K < 1/2, on the contrary, is a line of unstable fixed points: the smallest V will 
grow upon renormalization if K < 1/2, the system will "go to strong coupling" and develop a gap in the excitation 
spectrum. 

At the XY point, vs = vp = J and K — 1. For small values of the couplings Jz, J2, ^3 the exponent K is larger 
than 1/2 and the umklapp term is irrelevant, V —^ 0. This region corresponds to a spin liquid. The large distance 
behaviour of the correlation functions is characterized by power laws with exponents related to K. For instance, 
spin-spin correlations like (SqS^) are related (recall that = — 1/2) to density-density correlations of the spinless 
fermions (n(O)n(x)). The density operator has a continuum limit expression of the type 

n{x) - (p+ix) + p^{x)) + [^l{x)yj^{x)+H.c.] 

= -^V^{x) + — sin [V^^{x) + 2kFx] . (38) 

Using the fact that correlation functions of the bosonic field are simple to calculate for the Gaussian model {V — 0), 
i.e., 

1 

G{x) ^ {<i>{x)m - $2(0))v=o - ^ In ^^-^ , (39) 



^i7<I>(2;) -^7*(0)^ 



V=0 



7V(47r) 



(40) 



it is simple to show that 



(n{0)n{x)) 



K 



2tt^x^ 



A 



cos 2kpx 



r.2K 



(41) 



A being a non- universal constant. 

Increasing the values of the couplings, the spin liquid phase becomes unstable, at K — 1/2, against two different 
gapped phases, depending on the sign of the umklapp term V. For V > (large values of J^) the strong coupling 
fixed point is characterized by a field $ which is frozen at a value such that cos VIStt^ — —1, i.e., ^/Att^ = Tr/2. It is 
then clear that density-density correlations acquire LRO, since from Eq. 38 we get 



{n{0)n{x)) ^ sin {tt/2 + 2kFx) = cos (2fci?x) 

signalling a site-centered 2kF charge density wave (CDW). In the spin language this corresponds to a Neel phase with 
titi LRO. For < 0, on the contrary, the strong coupHng fixed point is characterized by a field $ which is frozen at 
the value $ = (or 2ti). To guess what correlation functions acquire LRO, notice that the canonical transformation 
il)p e~*P'^/^'0p changes the sign of the umklapp term in Eq. |2^ (in boson language this corresponds to \/4^^ 
^/Att^ + tt/2). Knowing that the 2A;i^--component of the density operator acquires LRO in the Nccl phase {V > 0), we 
immediately conclude that the operator having LRO for V < reads 



i[ijUx)ilj^{x) -H.c] 



1 



cos [v47r$(a;) -I- 2kFx] 



(42) 



An operator whose continuum limit 2/cF-component is given by Eq. |42 is readily found to be the bond kinetic energy 
(c|ci-f 1 + c|_j_j^Ci). The strong coupling phase obtained for < is therefore a bond-centered charge density wave, to 
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be contrasted to the site-centered CDW obtained for > 0. In spin language, this bond-centered CDW is a spin 
dimer phase. 

Separating the Neel {V > 0) from the dimer phase (V < 0) is the hne of unstable fixed points {V — with K < 1/2) 
mentioned above. Along this line ('PM' in Fig. the system behaves as an effective Luttinger (or Gaussian) model 
with 1/8 < K < 1/2. If i^T < 1/8, cosine terms of the type V cos - formally coming from higher order 

umklapp processes involving four-particle scattering - become relevant and open up a gap. The nature of the transition 
line changes to first order. Correlation functions behave as power laws along the line PM. Density-density correlations, 
for instance, are still given by Eq. All the critical exponents along the V = 0-line are known in terms of K. The 
correlation function exponent follows directly from Eq. i.e., t] = 2K. The gap between the ground state and the 
first excited state goes like |^ 

A = i - |y|i/(2-4if) ^ 

implying a correlation- length exponent v = 1/(2 — AK). The order parameter exponent is given hy P ~ vK. The 
specific heat exponent is a = 2 — 2i/ = (2 — 8K)/{2 — 4K). 



V. THE SPIN DIMER PHASE AND ITS SURFACE INTERPRETATION 



In the spin dimer phase, ordinary spin-spin correlations decay exponentially to zero, but four-spin correlation 
functions of the type {{Si ■ Si+i){Sj ■ Sj+i)) acquire LRO. More specifically, everywhere inside the dimer phase 
in Figs. HH, one has: 

5f " = ( (^o'^i^) (5|5|+i) ) « A + B{-iy j^oo. (43) 

This is illustrated in Fig. ^ where we show the size dependence of various static structure factors at the point 
{Jz = 3J, J2 = 2.4J, J3 = 0). These values are obtained from exact diagonalizations of chains up to 28 sites. The solid 
squares represent the dimer static structure factor at q = n 

S"i™(q = n)^J2 e^^'-Sf , (44) 
j 

whereas the open squares and the stars represent, respectively, the it/2 and tt component of the ordinary spin-spin 
structure factor. Clearly, S'^^"^{q — n) diverges linearly with the length of the chain (see inset of Fig. ||), implying 
long-ranged oscillations of the corresponding correlation function, whereas the usual spin-spin structure factor is finite. 

To illustrate in more detail some of the physics of this disordered spin state, and its translation into the surface 
language, we consider a representative dimer phase point. As it happens, there is a special point in the phase diagram 
of the J — J2-model (Eq. 22| with Jz = J — 2J2), where the twofold degenerate exact ground state is exactly known, 
p5| and extremely simple: it is just a product of spin singlets. Explicitly, for any finite (even) size N the two ground 
states, which turn into one another by translation of a lattice spacing, are 

1*1) = |12)|34) •••|iV- lA^) 

|*2> = |23)|45) •• - lA^ - 2A^- l)|Afl) . (45) 

Here \ij) = | ti ^ i^)/V2 denotes a singlet between sites i and j. Eq. ^ is just the spinless fermion translation of 
^'2. Some of the properties of these prototype dimer states which we are going to illustrate have been first discussed, 
in connection to the DF phase problem, in Ref. [ p^ . Obvious properties of l^^i) are, for instance, that spin-spin 
correlations are extremely short ranged, 

(*i|5||«'i)=0 Vj 

=0 \t-j\>l, (46) 

and that translational invariance is spontaneously broken, 

(*i|52Vi^ljl*i) = -1/4 

(*i|^2,52Vil*i) =0- (47) 
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In spite of this order parameter, such states are clearly spin- disordered. (1 x 1) order for the surface, for instance, 
translates into Neel LRO for the spin chain (see Fig. ||), whereas a dimer state has only short range spin-spin 
correlations. To see why they describe a flat surface, consider expanding the product of singlets in Eq. for l^*!), 
say. One obtains the sum of 2(^/2) spin configurations, one of which will be of the typical form 

(Ti)(Ti) (iT)(iT)(iT) (Ti)(Ti)-- - • (48) 

Here we have taken the (tJ,) part of the singlet for the first two pair of sites, the (J,t) part of the singlet for the next 
three pairs of sites, and so on. A down (2 x 1) step (i.e., a pair of neighboring down spins, see Fig. ||) is obtained 
each time a (J,|) pair follows immediately after a (ti) one, and, vice-versa, an up (2 x 1) step (a pair of neighboring 
up spins) results from a ("fl) pair following a (J,|) one. In between steps, there are regions with Neel type of order 
(unreconstructed regions in the surface language). Clearly, there is no way of having two up steps (or two down steps) 
following each other: a step up is followed necessarily by a step down and vice-versa. The surface is therefore flat. 

In the dimer phase there are characteristic correlations between steps that are worth stressing. An up (down) (2 x 1) 
step ending at site j is "measured" (see previous discussion) by the spin operator 

Step± = (5|_i±l/2)(5|±l/2) , 

counting, respectively, IT (Step^) and || (Step^) combinations at sites {j — l,j). An operator counting a step, 
irrespective of its being up or down, is given by 

Stepj- = Step+ + StcpT = 2{S^_^S^ + 1/4) . (49) 

One can easily work out correlation functions for such step operators in the representative dimers states. For odd j, 
for instance, one finds: 

r i if n = 
(^'i|Step+Step+^„|5'i) = \ je if n > 2 and even , 

[ otherwise 



i if n = 2 

(^'i|Step^StepJ_^„|*i) = "( ^ if n > 2 and even 

otherwise 



( i ifn = 

(^'i|StepjStepj+„|*i) = < I ifn>2andeven . (50) 

[ otherwise 

Similar results apply to ^2 for the case of even j. It is interesting to see how closely a point inside the dimer phase of 
Fig. H resembles such an ideal scenario. Fig. show step-step correlations (StcpjStepj^„) obtained from exact 
diagonalization of a chain of 28 sites, for a point inside (a) the dimer phase {J^ — 3.0J,J2 = 2.4J,J3 = 0), and (b) 
the Neel phase (J^ = 3.0J,J2 = J3 = 0). In the Neel phase, IT a-nd || steps are bound in pairs, and the correlation 
function decays exponentially to the square of the step-concentration, shown by a dashed line in Fig. |o| (b). The 
relevant defect is therefore the domain wall denoted by £2x1 in Fig. |^. In the dimer phase, on the contrary. If and 
a steps are unbound, and free to move in a fluid-like manner, but their correlation function displays long-ranged 
oscillations with period tt. In other words, the fluid of (roughly) alternating up and down steps has the feature that 
steps prefer to stay at an even distance from each other. In the neighboring (2 x 1) phase, this fluid of 2 x 1 steps 
solidifies into an ordered structure of the type ttii- We stress the fact that the oscillations displayed in Fig. |l^ (a) 
are not due to (2 x 1) order; the point considered is, as demonstrated in Fig. ^, disordered. 

Step-step correlations of the type shown are simple manifestations of the spontaneous breaking of translational 
invariance. Similar (and related) effects can be seen in other properties of the disordered surface. Suppose we want 
to count, in the surface terminology, the difference in the number of white and black local maxima in the surface. 
We restrict first our considerations to sites which are local maxima when considered in the x-direction only. In the 
spin language, a local "maximum" at site j occurs whenever the site j — 1 has spin "f and the site j has spin |. An 
operator which "counts" the maximum at j is therefore {Sj_i + l/2)(l/2 — S^). The difference between white (even 
j) and black (odd j) maxima is therefore measured by the order parameter 
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Pj;^''^ = {2/N)J2e^-'{S]_, + l/2)(l/2- . 
j 

^BW^^ is odd under translation. Its value is 1 on the Neel state | titi • • ^^'^ ~1 the other Neel state | J,tiT ' ' ■)■ 
Quite generally, it is different from zero in the whole Neel phase of the spin phase diagram. Consider now the value 

of P^^^'^^ on the dimer state Using the elementary results in Eq. p6[^ , we arrive at 

{^i\Pif^\^i) - -|^E'^"-''(*il^l-i^|l*i) = i • (51) 

i 

Similarly, (^2|-Pbh/"^I^2) = —1/4. Therefore, the implication of the dimer scenario, with its spontaneous breaking 
of translational symmetry, is that, on the disordered flat surface, one of the two sublattices tends to dominate in the 
local maxima. 

One can check this prediction by Monte Carlo simulations of the original classical models. In the next section we 
will present the results of our simulations for the and the K4-iiiode\. The results strongly support the dimer phase 
scenario. 



VI. MONTE CARLO RESULTS AND DISCUSSION 



We have performed classical Monte Carlo simulations of the A'3- and iir4-model in the DF phase. We have measured, 
to start with, standard quantities like the square mean width of the surface, Sh^, 

= (^E(^r-V)'), (52) 
r,r' 

the (1 X 1) order parameter, Pixi, and the (2 x 1) reconstruction one, P2XI, 

^ r ^ reW 

^xi^^^E'^re''''""^^). (53) 



P2 



Here Nc is the number of cells in each sublattice (i.e., 2Nc is the number of atoms), and G = (27r/aa,)x. The square 
mean width Sh^ diverges logarithmically in the rough phase as the size of the sample L 00 

dh^ w K{T)\ogL , 

with a coefficient K(T) larger than a (universal) minimum value KiTfs) — I/tt^ attained at the roughening temper- 
ature. Pixi is different from zero only in the unreconstructed region of the phase diagram and goes to zero at the 
preroughening line. P2XI is different from zero in the reconstructed region of the phase diagram and goes to zero at 
the "Ising" line. Clearly, the DF phase has Pixi = 0, P2XI = 0, and 5h^ < 00. On the basis of the spin mapping and 
of the dimer phase scenario we expect, however, that some form of order will be present: one should be able to tell 
which of the two sublattices (W or B) prevails in the top layer. A way of testing this is to define the "local peak" 
operator 



Or 



1 4 

-n[AV. + l] (54) 



16 - 1 



where Ah^^i = h^ — /ir+bi and with i — I, ■ ■ ■ ,4 are the vectors connecting a chosen site to its four nearest neighbors 
(belonging to the opposite sublattice). Or takes the value 1 for the atoms lying above all their neighbors, and zero 
otherwise. Summing over all the sites with a phase factor 1 for the W sites and —1 for the B ones, we get a quantity 
measuring which sublattice prevails in the top layer. 



l^E^^^'-Or)- (55) 
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As defined, Pbw is normalized to 1 on the unreconstructed ground states, and to 1/2 on the reconstructed (2 x 1) 
ground states. |37| ] Our expectation is that Pbw is different from zero in the disordered flat phase, and vanishes in 



the rough region and on the preroughening line. (See Fig. 15.) 



A classical grand-canonical single-move Monte Carlo code has been set up and used for lattices of linear size 
L = Nx — Ny up to 100. Starting from a disordered surface, we randomly add or remove particles, making sure that 
the BCSOS constraint is fulfilled at each step, and accept moves according to the standard Metropolis algorithm. 
The configurations resulting from consecutive sweeps of the lattice (2L^ attempted moves) are quite correlated, so 
that independent values for the various averages are obtained as a result of a sufficiently large number of Monte Carlo 
sweeps. It is on the basis of such "independent measurements" that statistical errors are estimated. Typically 20 to 
50 such measurements are performed, each of which consists of 10^ — 10^ sweeps, after a suitable equilibration of the 
system. 

For the Ks-model, we used the parameters of Mazzeo et ai, roughly chosen to fit the glue model results of Ercolessi 
et al. [Q for gold: K^xji^iy = —0.51, K'i/K2y — 0.22 (i.e., /C = —2.3). An Ising type deconstruction transition 
has been reported to take place at Td k, 2.90K2y, while a Kosterlitz-Thouless roughening transition has been found 
at Tr « 3.09i^2j,. flit] We have performed a careful finite-size scaling analysis of the different order parameters at 
the intermediate temperature T = 3.0K2y. The surface is still flat at this temperature, as demonstrated in Fig. [Tl| , 
showing that Sh^ versus logL stays definitely below the universal critical slope K{Tb) — I/tt^, which implies that 
SK^ will eventually saturate to a constant as L — > oo. Fig. |l2| (a) shows the results obtained for P2XI (solid circles), 
and Pbw (diamonds). The squares denote a further order parameter used by Mazzeo et al. ||lT[| , 



P, 



(2x1) _ 
\W ~ 



1^ 



(56) 



LrSW res 

where the classical "spin" variables S^. are defined in terms of the nearest neighbor height differences as 



S-r^VA/lr,,. (57) 



i=l 



P2XI vanishes as L ^, see inset of Fig. ^ (a), confirming that T = 3.0K2y is above the deconstruction temperature 
To, in agreement with Ref. [0. Both Pbw and Pbw^'' decrease, instead, much slower than L~^. Fig. |2| (b) shows 
the logarithm of Pbw versus logL. The data for small sizes {L up to 48) can be fit with a power law L^" "^^. For 
larger values of L, a crossover is seen to what is most probably an exponential convergence to a non zero limit for 
Pbw- In other words, systems up to L = 48 are still smaller than the actual value of a correlation length S^bWi so 
that a fictitious power law behavior is initially seen. A similar behavior is also found for P^^/ . 

The corresponding results for a point inside the DF phase of the iir4-model {Ki/K2y = 0.1, K2xl^2y — —0.056, 
and T / K2y = 2.3) are shown in Fig. ^ Entirely similar comments apply to this case. 

The point in the Ki-model phase diagram to which Fig. |l^ refers to, is in fact located close to the preroughening 
line. A typical snapshot of the way this disordered flat surface looks like at this temperature is shown in Fig. |lj. 
Strictly speaking we are in a parameter region where the classical ground state is (2 x 1) MR reconstructed and the 
most energetically favored defects are Ising walls (see table 1). However, Ising walls in their ideal form (see Fig. |^) 
are almost totally absent. What one finds, instead, are extended walls of the Ising type with a width of arbitrary 
length. These are nothing but large (1 x 1) unreconstructed regions lying between two opposite (2 x 1) steps. Such 
(2 x 1) steps, which are the very building blocks of a MR structure, are now free to move in a fluid-like manner with 
the only constraint that an up-step is followed by a down-step. Occasionally, sequences of up-down (2 x 1) steps gain 
positional order by "solidifying" in (2 x 1) MR regions which are, however, always of the same "color" (more precisely, 
black, for the phase illustrated in Fig. Overall, the surface seems to have as many black regions as white ones: 
the Pixi order parameter, which counts precisely the relative abundance of W and B (1 x 1) elementary cells, is small, 
and goes to zero in the thermodynamic limit. Correlations between steps, however, or, in more elementary terms, the 
fact that every (2 x 1) step always ends into a B top atom, result in the above mentioned feature of the absence of 
white MR regions, and are such that Pbw turns out to be different from zero, albeit small. Altogether Fig. |lj is a 
nice illustration of how a dimer disordered flat phase should look like. 

These features should be of some relevance in the context of surface scattering experiments. We discuss here the case 
of He scattering. In the kinematical approximation, and within a SOS framework, the intensity of the specular peak 
(parallel momentum transfer Q « 0) with perpendicular momentum transfer in the so-called anti-phase configuration 
is given by 

|(Ee^"''-ar)|'<5Q,o + iVsitcsfcsrx(Q) , (58) 
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where is an appropriate "shadowing factor" which takes into account the physical requirement that surface peaks 
scatter more than valleys. |]39| , |ll| The first term is a (Bragg) coherent contribution, proportional to the square of the 
order parameter and of the number of sites. The second contribution, due to incoherent terms, is proportional to the 
susceptibility of the order parameter and to the number of sites. For our BCSOS-type of model, in which is even 
in the W sublattice, and odd in the B sublatticc, one immediately concludes that e*'^''"' = e'*^ "" for any allowed height 
configuration. The coherent part of the specular anti-phase peak /'^°'*(Q = 0,qz = tt/o^) would therefore be identically 
zero if all the surface atoms were to scatter in the same way (ar = 1 for all r). In the opposite assumption that only 
the local peaks scatter efficiently {ar = 1 if r is a local peak, ar — otherwise), we obtain that /™''(Q — Q,qz = tt/a^) 
is exactly proportional to the square of the Pbw order parameter p2| 

r°'*(Q = 0,q,= n/a,) cx iV^^jPswl' ■ (59) 

Quite generally, for a reasonably large class of choices of shadowing factors a^, the breaking of translational invariance 
should guarantee that I'^°'^{Q = 0,qz — Tr/a^) is different from zero (albeit possibly small) in the DF phase considered 
here. (More precisely, this is so for all the shadowing factors which can be written in terms of local operators of the 
hr variables, whose correlation function is long-ranged in the DF phase.) 

Experimentally, a dimer-type of disordered flat phase would manifest itself with a rapid fall of the anti-phase 
scattering as the critical temperature is approached, followed by an intermediate temperature region, where the surface 
is in the disordered flat phase, in which a small coherent anti-phase scattering intensity survives. This situation is 
sketched in Fig. By normalizing the scattering intensity to its low-temperature value, a dip at the critical 
temperature should be observable even if one considers the total scattering intensity /(Q — 0,qz = ir/az), which 
includes the incoherent contributions. (Strictly speaking, these contributions are proportional to the susceptibility 
which diverges at the critical temperature as where L is the size of the system. Since = NsUcs, the 

incoherent contributions will never win over the coherent part (cx .^g^itcs) overall dip should be observable 

in the normalized scattering intensity at the critical temperature.) Clearly, an important requirement for the dimer 
scenario, which one should test experimentally, is that the dominant defects proliferating on the disordered flat surface 
are indeed monoatomic, or (2x1), steps. The correlations of such monoatomic steps are, at least in principle, also 
accessible by direct imaging techniques, such as fast STM. |Q 

We mention here, before ending the section, a particularly simple choice of shadowing factors, proposed in Ref. [ p9| , 
which does not involve long-ranged operators: 

where is the number of neighbors of the atom in r which are found at a level higher than the atom itself. This 
expression linearly interpolates between a = 2 (local maximum) and a = (local minimum), and can can be recast 
in the form 



I 4 

1=1 



where are the vectors connecting site r to the four neighboring sites. Indeed, by exploiting this linearity, it is very 
simple to show that such a choice of ar leads to a I'^°^{Cl = 0, = n/az) which is proportional to the square of the 



(1x1) order parameter Pixi (see Eq. 53), 



/-^Q = 0,g, = ^/a,) « |(^e'°-"-/ir)P = AfflPixiP , (60) 

r 

and therefore vanishes at and beyond the preroughening line. Therefore, experimental scattering geometries should 



be chosen so as to emphasize peak-atom scattering, if the non monotonic behavior of Fig. 15, typical of Pbw-, is to 
be detected. 



VII. CONCLUSIONS 



The motivation for the present work was a deeper understanding, based on well-defined hamiltonians, of the nature 
of the disordered flat phase (or phases) occurring in simple lattice models of fee (110) surfaces. In particular, for 
reconstructed surfaces, in the spirit of the distinction proposed in Ref. |l^ between a DEF (Ising wall dominated) 
phase as opposed to a DOF (step dominated) phase, we wished to clarify which of the two scenarios was at play in 
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simple BCSOS-type models. The outcome of our study is that neither of those simple prototypes apphes, strictly 
speaking, to the description of the disordered fiat phase we find, which is, on the contrary, closely related to dimer 
phase of one dimensional quantum spin-1/2 systems. 

The phase diagram in Fig. ||, very similar to the one discussed in Ref. , shows many features that we believe to 
be quite robust. First, a transition between the (2x1) MR reconstructed phase and the DF phase with exp onent s 
which appear to be very close to Ising. (Although the actual nature of the line is an open issue, see end of Sect. IV A .) 
Second, the transition line between the unreconstructed and the DF phase (preroughening) has variable exponents, 
as was predicted. Third, the disordered phase has a non-trivial order parameter Pbw- It is quite remarkable that 
both the microscopic models discussed here and the cell model of Ref. ^ point in the same direction, to a disordered 
flat phase which has a non- vanishing order parameter of the type of Pbw- We recognize that such a feature is also 
present in the phase diagram of Ref. |Q . 

The obvious open question is whether the disordered phase discussed above is the only one possible. In other words, 
can we build microscopic models were defects other than (2 x l)-steps play a role and the resulting disordered flat 
phase (or phases) has qualitatively different features ? 

The discussion has to consider separately, at this stage, the case of semi-microscopic cell-type models [|j from that 
of fully microscopic surface models. In the former framework of a coarse-grained description of the system, as the four 
state clock-step model of den Nijs the stage is clear and the actors are there: walls and steps. Since white atoms 
stay on top in regions where the reconstruction variable is or tt, and black atoms do so in regions where is 7r/2 or 
37r/2, the Pbw order parameter has to be non-zero in the disordered flat phase of this model, which could therefore 
be called DEF. Indeed, the disordering transition resulting in a DEF phase is mostly promoted by walls, which involve 
a change of tt for 9 on either side of the defect. On the contrary, Pbw is expected to vanish in a hypothetical DOF 
phase, since in this case the relevant defects are steps, which involve a change of ±7r/2 for 9. The four values of 
9 should appear with the same probability in such a disordered phase, and there is no way of telling which "color" 
prevails in the top atoms. 

The question of possibly finding a DOF phase in the sense of Ref. in a model of the clock-step type deserves, 
however, a few comments. Suppose that steps were indeed the most energetically favorable objects in the problem. 
Eg <g; Eu], and imagine desiring a stable DOF phase, i.e., preventing the appearance of steps from making the surface 
immediately rough. The natural way of doing this is to assign vertex energies to the crossing of steps, in such a way as 
to disfavor the crossing of parallel steps with respect to antiparallel ones, [ p^ This is indeed the standard mechanism 
by which a DOF is stabilized in the context of RSOS models for simple cubic (100) surfaces. Such 6- vertex energies 
have been neglected by den Nijs in deriving the zero-chirality limit phase diagram for the clock-step model. |^ It is 
therefore an appealing suggestion, deserving further study, that their proper inclusion might open up the possibility 
of a genuine DOF phase in the model. 

Microscopic SOS models are in many ways attractive, at first sight, as far as stabilizing a DOF is concerned: they 
automatically tend to disfavor crossing of parallel steps which involve large height differences. Moreover, tuning 
the model parameters offers, in principle, the possibility of making steps or walls more favorable, at least as far as 
their T = energy is concerned. Things are however not so straightforward in practice. Consider, as a remarkable 
counter-example, the case of the i4r4-model. When — 1 < K^xjK^ < 0, the ground state is (2 x 1) MR reconstructed, 
and simple Ising wall defects are energetically more favorable with respect to all kinds of steps (see table 1). A value 
of K2X/K4 = —0.56, which we considered in one of the simulations, would have seemed therefore a quite promising 
candidate for a DEF phase. What we end up with is, instead, a situation quite well represented by the snapshot in 
Fig. |l^. The state of the system looks as predicted for a dimer spin state. We clearly see that there are large regions 
in which the surface looks unreconstructed (with either the W or the B sublattice on the top layer), separated by 
(2 X 1) steps, forming a fluid with up-down order but without positional order. It is worth stressing that the relevant 
objects in such a disordered phase - the (2 x 1) steps - are the most natural defects of the unreconstructed surface, 
and the very building blocks of the neighboring MR reconstructed surface (which can be seen as a solid of alternating 
(2x1) steps). The result of tailoring the T — defect energies in such a way as to promote an Ising- wall dominated 
DEF phase, ends up with an amusing realization of a dimer state instead. 

In conclusion, we believe that a dimer phase type of disordered flat phase is a natural candidate in systems with 
a BCSOS-type of symmetry like the fcc(llG) surfaces considered in this work. Experimental signatures of such a 
scenario would be the the detection of a rapid fall of the anti-phase scattering intensity as the critical temperature is 
approached, followed by an intermediate temperature region (before roughening) in which the dominant defects are 
monoatomic steps, and where a small coherent anti-phase scattering intensity survives. 
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FIGURE CAPTIONS 



FIG. 1. Schematic top view of the fee (110) surface. The two sublattices, W and B, are denoted by open and solid circles. In 
the ideal unreconstructed (110) surface, one of the two sublattices lies at a distance = ay/2 above the other. The couplings 
considered in the model are indicated. Lattice basis vectors are also shown. The dashed zig-zag lines represent successive row 
configurations ( "time-slices" ) used in the spin-chain mapping. 

FIG. 2. Relevant extended defects (steps and walls) of a (1 x 1) and of a (2 x 1) reconstructed surface. The ground state 
energies of these defects are given in table 1. £2x1 is the (2 x 1) (or monoatomic) step discussed in section ^ (2 x l)-steps 
prohferate in the DF phase, mantaining up-down long-range order. £2x1 is a bound pair of (2 x l)-steps, the relevant defect of 
an unreconstructed surface, ecs and eAS are clockwise (or (3 x 1)) and anticlockwise (or (1 x 1)) steps, eiamg and e/si„g are 
two possible types of domain walls. 

FIG. 3. Schematic height profiles of the two ground states (Ul=White, U2=Black) of the unreconstructed surface, and of 
the four ground states of the (2 x 1) missing-row surface. The reconstruction variable 9 is indicated. The spin representation of 
each state, using Eq. ^ is explicitly given. Notice that the two unreconstructed ground states correspond to the two possible 
Neel states of a spin-1/2 chain. 

FIG. 4. The exact roughening temperature of the anisotropic BCSOS model (solid line) and the result obtained by making 
use of the mapping onto the XXZ Heisenberg chain (dashed line), as a function of the anisotropy ratio KixIKiy. The inset 
shows the relative discrepancy between the two results. 

FIG. 5. Ground state phase diagram of the Heisenberg chain with second neighbor SlSl^^ coupling. Ground state 
degeneracies are given in square brackets, and the translation of the different phases in the surface language is explicitly 
indicated. The (1 x 1)-DF line starting at the point P is continuous, with a variable exponent. 

FIG. 6. Qualitative ground state phase diagram for the Heisenberg chain with second- and third- neighbor couphngs, 
for J3 = QJ2 with < a < 1/2. Ground state degeneracies are given in square brackets. The line labelled 'PM' has a variable 
exponent. 



FIG. 7. Phase diagram for the Jfs-model, as obtained from the quantum spin-chain mapping, for Ki/K2y ~ 0.025. 



FIG. 8. Qualitative phase diagram for the modified BCSOS-model in Eq. yj for fixed values of the couplings K3 and K4. 

FIG. 9. Finite size behavior of different structure factors at the point Jz/J = 3.0 and J2/J — 2.4 in the dimer phase: 
S"^'™(7r) (solid squares), S^^{it) (crosses), and S^^{it/2) (open squares). The inset shows a logarithmic plot of the dimer 
structure factor, together with a dashed fine of slope 1, for comparison. 

FIG. 10. (2 x 1) step-step correlations, see Eq. ^ for a chain of 28 sites, (a) at the point (J^ = 3.0J,J2 = 2.4J,J3 = 0) 
inside the dimer phase, and (b) at (J^ — 3.0 J, J2 — J3 — 0) inside the Neel phase. For the ideal dimer state 'Jfi, see Eq. ^ the 
correlation function would oscillate between the values and 1/4 (the dashed line in (a)). Inside the Neel phase the correlation 
function decays exponentially to the square fo the step-density, denoted by the dashed line in (b). 

FIG. 11. Finite size scaling of the height fluctuations, Eq. ^ for the K3 model at K2x/K2y ~ —0.51, Ks,/K2y = 0.22, 
T/K2y = 3. A line with the critical slope K{Tr) = I/tt^ is also shown, indicating that the surface is smooth at this point. 
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FIG. 12. (a) Finite size scaling of P2XI (the reconstruction order parameter, Eq. ^ full circles), Pbw (Eq. open 
diamonds), and Pg^^'' (Eq. open squares) for the K3 model at the same point considered in Fig. |l^. The inset shows that 
-P2XI vanishes as the inverse of the linear size L of the lattice. The surface is thus deconstructed, (b) Log-log plot of the finite 
size behavior of Pbw, showing the saturation to a non-zero value for L — > 00. The surface is in a disordered flat state. 

FIG. 13. Same as in Fig. ^| for a point inside the disordered flat phase of the i('4-model {K2x/K2y — —0.056, K4,/ K2y ~ 0.1, 
T/K2y = 2.3). 

FIG. 14. Snapshot of a surface configuration as generated by the Monte Carlo simulation for the 7<'4-model at the same 
point considered in Fig. inside the (dimer) disordered flat phase. 

FIG. 15. Sketch of the expected behaviour of the Pbw order parameter (Eq. ^), proportional to the antiphase scattering 
intensity, as a funtion of temperature when the preroughening line is crossed. 
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TABLE 





J^s-model 


_ft'4-model 


£2x1 


4/s:2x 


4K2:c + 8K4 


^2X1 




8K2:c + 8K4 


ecs 


2K2a: + 8K3 


2K2^ + 8X4 


^AS 




-2K2:, + 8Ki 


Rising 




-4K2, + 8K4 




4K2^ + I6K3 


4.K2^ + I6K4 



TABLE I. Ground state energy of the defects shown in Fig. 2 for the K3- and i(r4-model. 
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(2x1) MISSING-ROW GROUND STATES 
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